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Abstract 

^ . We present the one-loop QCD amplitudes for two external massless quarks and three external 

gluons (qqggg) ■ This completes the set of one- loop amplitudes needed for the next-to-leading-order 
corrections to three-jet production at hadron colliders. We also discuss how to use group theory 
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one- loop two-quark n-gluon amplitudes. We use collinear limits to provide a stringent check on the 
amplitudes. 
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1. Introduction 



Jet physics at hadron colliders allows one to confront the theoretical predictions of QCD with 
experimental results and thereby probe for new physics at the highest possible energies. Yet precise 
comparisons between theory and experiment are hampered by the lack of calculations beyond the 
leading order of perturbation theory, for all but the simplest processes. In pure QCD, the next- 
to-leading-order corrections computed to date [1] have relied on the one-loop amplitudes for four 
external partons, first calculated by Ellis and Sexton [2]. More recently, we have calculated the 
one-loop amplitudes for five external gluons (ggggg) [3], and Kunszt, Signer, and Trocsanyi (KST) 
have calculated the amplitudes for four quarks and a gluon (qqqqg) [4]. In this paper we present 
the remaining one- loop five-parton amplitudes, for two (massless) quarks and three gluons (qqggg)- 
Combining these analytic results with the known six-parton tree amplitudes [5,6], one can now con- 
struct numerical programs for next-to-leading-order corrections to three-jet production at hadron 
colliders, and examine the structure of jets, for example dependence of cross-sections on the cone 
size, beyond the leading non-trivial order probed in next-to-leading order two-jet programs [1,7]. 
Computation of the ratio of three-jet to two-jet events at hadron colliders at next-to-leading order 
in a s would also make possible the measurement of a s in purely hadronic processes and at the 
largest energy scales available. 

Many methods developed in recent years can be used to simplify the computation of one- 
loop multi-parton amplitudes, including spinor helicity methods [8], color decomposition of am- 
plitudes [5,9,10], string-based techniques [11,12,13,14,3], supersymmetry Ward identities [15,16], 
supersymmetry-based decompositions [3,17,18], and perturbative unitarity [19,20,21]; all of these 
techniques have been used to obtain the amplitudes presented in this paper. 

We have found it useful to organize the calculation in terms of gauge-invariant, color-ordered 
building blocks, dubbed primitive amplitudes. We show in the next section that all of the kine- 
matic coefficients (partial amplitudes) appearing in the color decomposition of amplitudes with two 
quarks and (n — 2) gluons can be expressed as sums over permutations of gauge invariant primitive 
amplitudes. The analytic structure of a primitive amplitude is generally simpler than that of a 
partial amplitude; a primitive amplitude receives contributions only from diagrams with a fixed 
ordering of external legs, while the generic partial amplitude receives contributions from multiple 
orderings. Thus, fewer kinematic invariants appear in each primitive amplitude. Although this 
organization was motivated in part by string theory, our discussion is entirely field-theoretic. 
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We use supersymmetry to reduce the number of quantities to be calculated. QCD ampli- 
tudes may be decomposed in terms of supersymmetric and non-supersymmetric parts. Through 
use of supersymmetry Ward identities, the supersymmetric parts of amplitudes with two external 
quarks and three external gluons may be obtained directly from the previously calculated five-gluon 
amplitudes [3]. 

We have also made use of the cut-reconstruction method described in refs. [20,21]. If certain 
power-counting criteria are satisfied, amplitudes are entirely constructible from their cuts. Al- 
though QCD amplitudes are generally not cut-constructible, by taking linear combinations of QCD 
amplitudes with ones involving scalars and/or gluinos, the QCD amplitudes may be separated into 
cut-constructible and non-cut-constructible parts. We have used this unitarity-based technique to 
obtain the cut-constructible components of some of the primitive amplitudes for qqggg (those that 
enter into the subleading-in-color contributions to the virtual part of the cross-section). Here the 
cut-constructible components are formed by adding to the desired diagrams a new set of diagrams, 
which differs only in the replacement of virtual gluons in the loop by scalars. For a specific choice of 
the Yukawa coupling between the scalars and the quark line, the sum of gluon and scalar diagrams 
satisfies the power-counting criteria (see ref. [21]). We then calculate the scalar contributions di- 
rectly; they are not cut-constructible, but they are easier to calculate directly than the full gluon 
contributions. Finally we reassemble the desired gluon contributions. 

In order to ensure the correctness of the amplitudes, we have performed a number of checks. As 
the momenta of two external legs become collinear the amplitudes must factorize properly. We have 
verified this factorization for all amplitudes in all channels. This provides an extremely stringent 
constraint on the amplitudes. In fact, this constraint is sufficiently powerful that it has been used 
to construct ansatze for a number of amplitudes with fixed helicities but an arbitrary number of 
external legs [22,23,20], which were then proven correct by either recursive [24] or unitarity [20,21] 
techniques. (The recursive and unitarity techniques have also been used to construct a variety of 
other one-loop amplitudes with an arbitrary number of external legs [24,20,21].) 

We have performed additional checks on certain helicity amplitudes by computing all diagrams 
that enter into a supersymmetry Ward identity [15], and explicitly verifying the identity. Not only 
does this provide a check on amplitudes presented in this paper, but also on the supersymmetric 
combinations of the five-gluon amplitudes presented in ref. [3]. (A similar supersymmetry check 
using the five-gluon amplitudes has been carried out [25] for the qqqqg amplitudes reported in 
ref. [4].) As a final check, we have verified that the cuts in some amplitudes obtained by more 
direct diagrammatic means are consistent with unitarity. 

In section 2, we give the SU(N C ) color decomposition for amplitudes involving two external 
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quarks and n — 2 external gluons, as a sum of color factors multiplied by partial amplitudes. We 
also give a formula for the sum over colors of the interference between tree and one-loop qqggg 
amplitudes, in terms of partial amplitudes; this formula is required for the virtual part of the 
color-summed parton-level cross-section. The primitive amplitudes, which form the gauge-invariant 
building blocks for the amplitudes, are described in section 3. The precise relation of the primitive 
amplitudes to the partial amplitudes is given in section 4. In section 5 we give the main results of 
the paper, the primitive amplitudes for qqggg. Section 6 contains our conclusions. Four appendices 
contain technical details related to color algebra and collinear checks. Appendix I provides a 
derivation of the relation between primitive and partial amplitudes. Appendix II collects the one- 
loop four-point amplitudes [16,11] that appear in collinear limits of qqggg amplitudes, namely gggg 
and qqgg. Appendix III then illustrates the procedure for carrying out collinear checks, using these 
amplitudes and "splitting amplitudes" from ref. [20]. Finally, appendix IV shows how to use the 
two-quark (n — 2)-gluon primitive amplitudes to construct amplitudes where some of the gluons 
are replaced by photons. 



2. Color Decomposition for Two-Quark (n — 2)-Gluon Amplitudes 

In this section we describe a color decomposition of the one-loop two-quark (n — 2)-gluon 
amplitude qqg . . . g, in terms of group-theoretic factors (color structures) multiplied by kinematic 
functions called partial amplitudes. In the following sections, we shall give formulae for all of 
the partial amplitudes in terms of color-ordered, gauge-invariant building blocks called primitive 
amplitudes. A primitive amplitude is defined as the sum of all one-loop diagrams in which the n 
external legs have a fixed order around the loop (the color order) , with some additional restrictions 
to be described in the following section. 

For the qqg . . . g amplitudes, let particle 1 be an antiquark, transforming in the N c represen- 
tation of SU(N C ), with color index i±, and let particle 2 be a quark, transforming in the N c with 
index «2- Denote these particles by lq and 2 q in order to distinguish them from the remaining 
gluons, particles 3 to n, transforming in the adjoint representation with indices 0,3, . . . , a n . We also 
allow for rif Weyl fermions and n s complex scalars circulating in the loop, both in the (N c + N c ) 
representation {rif flavors of massless quarks and n s massless scalars). 

The general strategy for obtaining multi-parton color decompositions is to rewrite the SU(N C ) 
structure constants f abc in terms of the group generators in the fundamental representation T a , 
normalized so that Tr (T a T b ) = 5 ab , 

jabc = — 1 Tv([T a ,T b ]T c ) . (2.1) 
V2 
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Then one applies the SU(N C ) Fierz identity 

(X 1 T a X 2 ) (Y 1 T a Y 2 ) = (X 1 Y 2 ) (Y 1 X 2 ) - ^{X 1 X 2 ) (Yi Y 2 ) , (2.2) 

where Xi,Yi are strings of generator matrices T ai , in order to remove contracted color indices. 
In discussing amplitudes where all external particles are in the adjoint representation (such as 
amplitudes in supersymmetric QCD with no matter content), or trees consisting of only particles 
in the adjoint representation, one may (and should) replace the SU(N C ) Fierz identity with the 
corresponding U(N C ) identity, since it is simpler: the 'photonic' term decouples [5,9,10]. For the 
two-quark amplitudes under consideration, at tree level only one string survives and the color 
decomposition is [5] 

< ee (l 9 -,2 9 ,3,...,n) = g n ~ 2 £ (T a ^ . . . T a ^)£ A^°(l q , 2 q ; <r(3), . . . , a(n)) , (2.3) 

where S n - 2 is the permutation group on n — 2 elements, and A^ ee are the tree-level partial ampli- 
tudes. They are identical to the tree-level partial amplitudes for the process with gluinos replacing 
quarks, and are thereby related to the tree- level all-gluon amplitudes by supersymmetry Ward 
identities [15]. We adopt throughout the convention that all momenta are taken to be outgoing. 
Because the color indices have been stripped off from the partial amplitudes, there is no need to 
distinguish a quark leg from an anti-quark leg; charge conjugation relates the two choices. 
At one loop an additional trace may survive, and the color decomposition is 

n-1 

A n (l q ,2 q ,3,...,n) = g n J2 E Gr SV( 3 > ■ ■ ■ ,n)) A n;j (l q ,2 q ;a(S, . . . ,n)) , (2.4) 

where the color structures Gr^ are defined by 

Grg>(3,...,n) = N c (T°» . . .T-)^ , 
G^. q 2 \3-A,...,n) = 0, 
Grg ) (3,...,j + l;i + 2,...,n) = Tr(T° 3 . . . T aj+1 ) (T a ^ . . . T a ")£ , j = 3,...,n-2, 
Gr^l^S,...^) = Tr(r«3...r«") , 

(2.5) 

and S n -j = Zj-i is the subgroup of S n - 2 that leaves invariant. When the permutation a 

acts on a list of indices, it is to be applied to each index separately: <r(3, ... ,ri) = c(3), . . . , cr(n), 
etc. We refer to A n;1 as the leading-color partial amplitude, and to the A n; j >1 as subleading, 
because for large iV c , A n -i alone gives the leading contribution to the color-summed correction to 
the cross-section, obtained by interfering »4^ ree with A n . The explicit N c in the definition of the 
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leading color structure Grj^ — which is otherwise identical to the tree color structure — ensures 
that A n -i is 0(1) for large N c . 

It is useful to recall the analogous color decomposition for n external particles in the adjoint 
representation, in particular for the pure super- Yang-Mills amplitude for two external gluinos and 
n — 2 gluons (similar to the decomposition for n-gluon amplitudes [10]), 



|n/2J+l 

Al VSY (l § ,2 § ,3,...,n) = g n E Gr nii ( < 7(l 5 ,2 5 ,3,...,n))^S SY ( < 7(lg,2 5 ,3,...,n)) ) 

j =1 ves n /s n;j 



(2.6) 



where 



Gr„.i(l,2,...,n) = N c Ti(T ai T a2 . . .T a ") , 

(2.7) 

Gr n;j (l,...,j-l;j,...,n) = Tr(T ai ...T a ^) Ti(T aj . . . T a ") , j = 2, . . . , [n/2\ + 1, 



[x\ is the largest integer less than or equal to x, g stands for the gluino legs, and S n -j is the 
subset of S n that leaves Gi n -j invariant. The similar structure of the gluino amplitudes and quark 
amplitudes will help in understanding how to simplify the organization of the latter, particularly 
the subleading-color contributions, A n; j > i. 

The partial amplitudes ^4^ SY for ggggg can be obtained from five-gluon results [3] and su- 
persymmetry Ward identities [15]; these can in turn be used to reduce the work required in the 
quark case. Throughout this paper, we consider only supersymmetric amplitudes with no matter 
content. The use of supersymmetry in loop amplitudes implies a calculation using a super symmetry- 
preserving regulator, such as dimensional reduction [26] or the four-dimensional helicity scheme [11] 
(which are very closely related at one loop). To obtain results for the quark amplitudes in other 
schemes, such as the conventional dimensional regularization method [27] (often called the 'MS' 
scheme in the literature), one must shift the partial amplitudes presented here by a quantity pro- 
portional to the tree amplitude; the constant of proportionality has been determined by Kunszt, 
Signer, and Trocsanyi [16] (see eqn. (5.5)). 

Unlike the all-external-gluon case, we can subdivide the two-gluino (n — 2)-gluon color struc- 
tures further, depending on the trace to which the fermion charge matrices belong. Define Gr* ;j - 
and Gr^ ;j - to be respectively Gr n;j - with fermion charge matrices in the same and different traces; 
furthermore, require for Gr* ;j - that the fermion charge matrices lie in the second trace, which 
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increases the number of such structures from \n/2\ + 1 to n — 1. We may then write 



n-l 



A^ SY (l § ,2 § ,3,...,n)=g n J2 E G< ; 2 5 , 3, . . . ,n)) ^ SY (a(l § , 2 5 , 3, . . . ,n)) 



j = l crGS s ■ 



L»/2J+1 

+ 5" E E G< i ( < 7(l g ,2 g ,3,...,n))^S SY ( < 7(l5,2 5 ,3,...,n)), 

(2.8) 

where • and S^ ;J - are the sets of permutations of n elements that act nontrivially on Gr* ;j - and 
Gr^ ;j - respectively, and preserve the assignment of fermion charge matrices to the two traces. In ap- 
pendix I, we will find this separation useful when relating adjoint representation partial amplitudes 
to fundamental representation partial amplitudes. 

Returning to the quark case, the explicit decomposition of the qqgg one-loop amplitude is 



.44(15,2,, 3, 4) = g 4 



N c {T a *T a *)£ A 4; i(l 9 -,2 9 ;3,4) + N c (T^T a %^ A 4il (l q , 2 q ; 4, 3) 

(2.9) 

+ Tr(T a3 T a4 )«5 i2 11 A 4;3 (l 9 -,2 9 ;3,4) 



in agreement with the decomposition used by KST [16] (apart from notational differences such as 
the explicit factor of N c in Gr^- 1 and the ordering of external legs). 

The decomposition of the qqggg amplitude is [28] 



A(l g -,2 g ,3,4,5) = g 



N c ^(nor-wr-o)^ A 5;1 (l ? -,2 g ;a(3),(7(4),(7(5)) 
<res 3 

+ Tr (T a ^T a ^) (r-'W),* 1 A 5;3 (l q , 2 q ; a(3), <t(4); <t(5)) 



+ Tr(T a3 T a4 T a5 )C A 5;4 (l 9 -,2 9 ;3,4,5) + Tr(T a5 T a4 T a3 ) A 5;4 (l g -, 2,; 5, 4, 3) . 

(2.10) 

In the partial amplitude an additional semicolon separates the gluon sandwiched between the 
quark indices (the last gluon in A 5 - 3 ) from the other two gluons. 

The virtual part of the next-to-leading order correction to the parton-level cross-section is 
given by the sum over colors of the interference between the tree amplitude 4^ ree and the one- 
loop amplitude A n . Using the color decompositions (2.3) and (2.4), and the Fierz rules (2.2) it is 
straightforward to evaluate this color-sum in terms of partial amplitudes. Here we give the formula 



for the five-point case, qqggg, 



N 2 — 1 

2 Re [A iee *M = 2g 8 c 

colors 



(iV 2 -l) 2 A5 ; i(l 5 ,2, ;< 7(3), <t(4), <7(5)) 

- (iV c 2 - 1) (a 5;1 (Lj, 2,; <7(4), a(3), (7(5)) + A 5;1 (l 9 -, 2 g ; a(3), a(5), a(4)) 

- A 5 ; 3 (l,-, 2,; <7(4), a(5); a(3)) - A 5;3 (l q , 2 q ; a(3), a(4); a(5))) 
+ (iV c 2 + l) A 5;1 (l,-,2,;<7(5),<7(4),<7(3)) + (iV 2 - 2) A 5;4 (l g -, 2 g ; a(3), a(4), a(5)) 
+ A 5;1 (l ? -,2 9 ;a(4),(7(5),(7(3)) + A 5;1 (l 9 -, 2,; a(5), a(3), a(4)) 

- A 5 ; 3 (l 5 , 2,; a(5), a(3); a(4)) - 2 A 5;4 (l 9 -, 2 q ; a(5), a(4), a(3)) , 



(2.11) 



where the sum is over all permutations of the three gluons. 



3. Primitive Amplitudes 

In this section we introduce a set of gauge- invariant, color-ordered building blocks, which 
we call primitive amplitudes, that suffice to determine all the two-quark (n — 2)-gluon partial 
amplitudes A n .j. Explicit expressions for primitive amplitudes tend to be much more compact 
than those for the generic partial amplitude, because the legs are ordered. Only a subset of the 
kinematic invariants (^ ki) 2 appear as arguments of logarithms or dilogarithms in a color-ordered 
set of diagrams, namely those where all the momenta are adjacent with respect to the ordering, 
and this leads to simpler analytic structure for such objects. 

It is not obvious a priori that every partial amplitude can be expressed in terms of primitive 
amplitudes, because the generic partial amplitude receives contributions from diagrams with several 
different cyclic orderings of the external legs, and it does not receive contributions from certain 
classes of diagrams present in the primitive amplitudes. Nevertheless, it is possible to write every 
two-quark (n — 2)-gluon partial amplitude as a sum over permutations of primitive amplitudes; 
one can show (using string-inspired groupings of diagrams) that the unwanted diagrams cancel 
out in the sum. We perform the necessary SU(N C ) group theory manipulations in the double-line 
formalism [29]. The color decomposition for one-loop qqg . . . g amplitudes is analogous to that 
presented in ref. [10] for one-loop amplitudes with n external gluons, and the manipulations used 
to simplify subleading color contributions are similar to ones used in ref. [20]. We also employ 
supersymmetry [15,17,16], in order to reduce the number of primitive amplitudes that have to be 
calculated directly for qqggg. 
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To present the primitive amplitudes, we will find it convenient to use the language of color- 
ordered diagrams. These are basically Feynman diagrams from which color indices have been 
stripped, the relevant signs encoded in the topology of the graph. The vertices are no longer 
symmetric under exchanges of legs, and thus the external ordering of the legs becomes important: 
two diagrams which would be equivalent as Feynman diagrams are generally no longer so as color- 
ordered diagrams. The reader may find a description of tree-level color-ordered Feynman rules in 
refs. [6,17]. One obtains them by using the trace representation (2.1) of the structure constants f abc , 
as well as eqn. (2.2) or its U(N C ) counterpart for adjoint states. One then dresses the Feynman 
diagrams using the doubled color line notation for the adjoint-representation propagators, and 
single color lines for fundamental-representation ones [29]. The two terms in equation (2.1), and in 
equation (2.2), mean that one Feynman diagram can generate many different color-flow diagrams. 
By convention, we draw diagrams so that the ordering of legs is clockwise around the loop. 

In computing the coefficient of a particular color structure, or configuration of generator matri- 
ces T ai , one can then remove group theory factors from the vertices, arriving at the color-ordered 
rules shown in fig. 1, where the gluon Lorentz indices are [i,v,p,\ and outgoing momenta are 
k,p,q, and the fermion g is taken to be in the adjoint representation. The cyclic ordering of the 
legs in fig. 1 is important since there is a relative sign between the two orderings of the three point 
vertices; if one interchanges any two of the legs the vertex changes by a sign. This sign follows from 
the antisymmetry of the structure constants f abc . 




H V n 
WW = -i%f 



V2 



Figure 1. The color-ordered Feynman rules have antisymmetric three- vert ices. Straight 
lines represent fermions, and wavy lines gluons. 
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In various steps of our explicit calculations it is advantageous [13,17,18] to use a different 
gauge than the standard Feynman gauge used in fig. 1, such as background-field gauge [30] and 
Gervais-Neveu gauge [31]. In the various gauges the color-ordered three- and four-gluon vertices 
look different. However, partial amplitudes and primitive amplitudes are gauge invariant (we prove 
the invariance of the latter at the end of this section); therefore the formulae below, expressing 
partial amplitudes as sums of primitive amplitudes, which are derived using the rules in fig. 1, will 
hold in any gauge. 

If the fermion is in the fundamental representation, the same rules hold, but one of the three 
color lines at the vertex, the one flowing along side the fermion line, should be removed or 'stripped 
off'. This procedure of 'color-stripping' equates a contribution to an amplitude with external 
fundamental-representation fermions to a contribution to an amplitude with external adjoint- 
representation fermions. The latter differ from the former in a way we shall make precise in 
appendix I. 

In order to motivate the precise definition of primitive amplitudes for two quarks and (n — 
2) gluons, we first discuss amplitudes with external gluons only. Consider the set of diagrams 
contributing to the leading all-gluon partial amplitude A n -i(l, 2, . . . , n). The only diagrams that 
contribute are those with a single ordering 1 . . . n of legs around the loop, the ordering matching 
the associated color structure Tr(T ai T a2 • • • T a ™), as depicted in fig. 2. The color-ordered diagrams 
for this partial amplitude may be further distinguished by whether a gluon, a fermion, or a scalar 
circulates around the loop; each contribution is separately gauge- invariant because the coefficients 
of rif and n s in the full amplitude must be gauge-invariant. We denote them by A^\, with J the 
spin of the circulating particle, J = 1, |, 0. 




Figure 2. The external legs of color-ordered diagrams have a fixed clockwise ordering. 

The subleading n-gluon partial amplitudes A n -j > i (which are only present when an adjoint 
particle, not a fundamental, circulates in the loop) generically receive contributions from diagrams 
with many different cyclic orderings. However, it is possible [20] to express the A^)- as a sum 
over permutations of the A^ , as we review in appendix I. Therefore the partial amplitudes A^ 
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suffice to construct the full n-gluon amplitude. On the other hand, other than separating the 
contributions of internal particles of different spin, it is not possible to find gauge-invariant subsets 
of the diagrams that contribute to A^} v When all external legs are gluons, then, the are the 
irreducible gauge-invariant pieces of the amplitude, and serve as primitive amplitudes. 

If we were interested only in the two-gluino (n — 2)-gluon (supersymmetric) amplitude we 
could take the primitive amplitudes simply to be the leading-color partial amplitudes A n; i, just 
as in the n-gluon case, because the same subleading-from-leading permutation formula that holds 
in the n-gluon case applies to mixed gluino-gluon amplitudes as well. However, for the qqg . . . g 
amplitudes, one has to divide the sets of color-ordered diagrams into finer pieces before such an 
approach can succeed. The need for a finer division of the diagrams arises from the color flow 
which goes solely in one direction along a quark line, but in both directions along a gluino line. 
The gluino amplitudes A n .± can be decomposed further in a gauge-invariant way; the pieces of this 
decomposition are the primitive amplitudes, which are also the pieces out of which we may form 
the amplitude for external quarks. 

To ferret out these gauge-invariant subparts of the supersymmetric amplitude, we must trace 
the external fermion line through the diagram. Since the fermion line has an arrow, we can dis- 
tinguish one-loop diagrams according to which side of the loop the fermion line is on. We define a 
'left' class of diagrams to be those where, following the arrow, the fermion line passes to the left of 
the loop; the remaining diagrams are in the 'right' class. For example, in fig. 3a the fermion line 
enters the loop and turns left, passing to the left of the loop, so this is a 'left' diagram; fig. 3b is 
a 'right' diagram. In fig. 4a the external fermion line does not actually enter the loop, but still it 
passes to the left, so this is a 'left' diagram; fig. 4b is a 'right' diagram. In diagrams of the type 
shown in fig. 4 the particle circulating in the closed loop may be either a gluon, a fermion or a 
scalar; the left/right designation is applied in the same way. (The external fermion line is always 
used to make the distinction, not the fermion that might be in the closed loop.) 




Figure 3. In diagram (a) the fermion line (following the arrow) turns 'left' on entering the 
loop, and in diagram (b) it turns 'right'. 
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(a) (b) 

Figure 4. In diagram (a) the external fermion line passes to the 'left' of the loop, and in 
diagram (b) it passes to the 'right'. A gluon, fermion or scalar may circulate in the closed loop 
(hatched region). 

The 'left' and 'right' diagrams have to be separated from each other when the external fermions 
are in the fundamental representation because their color factors are different. (For gluinos in the 
adjoint representation the color factors are identical, and in fact the 'left' and 'right' diagrams 
must be added together to get supersymmetric partial amplitudes.) We shall show below that this 
division into 'left' diagrams and 'right' diagrams is gauge invariant; those 'left' and 'right' diagrams 
with a closed fermion or scalar loop form further gauge-invariant sets. We thus take the primitive 
amplitudes for qqg ... (7 to be 

^>[ J l(l g -,3,4,...,2 g ,...,n), 
R,[J] , ^=1.5.0, (3.1) 

A-n (Igj 3, 4, . . . , 2 q , . . . , 7l), 

corresponding to the sum of all diagrams with the indicated cyclic ordering of external legs, where 
the fermion line from lq to 2 q passes to the left (L) or to the right (R) of the loop, and where 
J = \ (J = 0) represents the subset of diagrams with a closed fermion loop (closed scalar loop). 
The normalization is such that two helicity states (Weyl fermions or complex scalars) circulate in 
the loop. Diagrams without closed fermion or scalar loops are assigned to J = 1; they may or may 
not contain a closed gluon loop, as the two types of diagrams mix under gauge transformations. 
We shall often suppress the superscript "[1]"; this creates no ambiguity. In the next section and in 
appendix I we show that all of the quark-gluon partial amplitudes A n; j can be obtained from these 
primitive amplitudes. 

Not all of the primitive amplitudes (3.1) are independent. By flipping over the set of diagrams 
where lq turns right, one obtains (up to a sign) the set of diagrams where lq turns left, with the 
cyclic ordering also reversed, 

^>[ J ](l 9 -,3,4,...,2 9 ,...,n-l,n) = {-l) n A^ J \l q , n, n - 1, . . . , 2„ . . . , 4, 3). (3.2) 

Also, the super- Yang-Mills partial amplitudes for two gluinos and n — 2 gluons ^ USY = ^4^ SY 
are given by the sum (with all cyclic orderings identical) 

4SUSY = .SUSY _ a L 1 a Ft > aL, [1/2] , aB.,[1/2] /o o\ 

J± n — i1 n;1 — J± n -f Ji n i- Ji n -f Ji n , iO.OJ 
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because the 'left' and 'right' diagrams have the same group-theory weight for an adjoint-representation 
fermion. In this equation supersymmetric cancellations occur between the 'left' and 'right' primitive 
amplitudes; in general, ^ USY is a much simpler quantity than either A% or A^. Equation (3.3) 
allows one to obtain one of the four terms on the right for free (say ^4^), given Ai^p^. 
Finally, the following fermion-loop contributions vanish, 

A^ 2 \l q ,2 q ,3,4,---,n) = A^(l q ,n,...,4,3,2 q ) = 0, 

(3.4) 

A^(l q ,3,2 q ,A,...,n) = A^(l q ,n,...,A,2 q ,3) = 0, 

and similarly for the scalar-loop contributions. The restriction to 'left' or 'right' diagrams combines 
with the ordering of the external legs to leave only tadpole and massless external bubble diagrams 
behind; but these are zero in dimensional regularization. 

We conclude this section by proving that the primitive amplitudes (3.1) are gauge-invariant.''' 
Suppose that the two external fermions (/i and $2) are cyclicly separated from each other by n a 
gluons going one way around the loop, and by gluons going the other way around (n a +rib = n—2), 
and set aside for a moment the closed fermion- or scalar-loop contributions. The sum A^ + A% 
is obviously gauge invariant because it equals the gauge-invariant, color-ordered partial amplitude 
^n!i SY ~~ ^n-i^ > i- e - the coefficient of the leading color structure 

N c Tr {T^T ai ■ ■ ■ T an *T f2 T bl ■ ■ ■ T bn i>) (3.5) 

in a theory where the fermion line is in the adjoint representation, and omitting the fermion loop 
contribution. (Different color structures are orthogonal in the large iV c limit, so their coefficients, 
when independent of iV c , must be individually gauge invariant [6].) 

To prove that A^ and A^ are each invariant independently, we consider a different gauge 
theory, with gauge group SU (N a ) x SU (Nb), where the n a gluons belong to SU(N a ), the rib gluons 
belong to SU(Nb), fermion 2 belongs to the representation (N a , iVj>), and anti-fermion 1 belongs to 
the representation (N a ,Nb). The coefficients of the color structures 

N a (T a "» • • • T ai )£ {T bl ■ ■ ■ T b ^ (3.6) 

and 

N b (T a "» • • • T ai ) I 1 (T bl ■ ■ ■ T bn » ) .f 1 (3.7) 

must be separately gauge invariant, since N a and Nb can be chosen independently. But the coeffi- 
cient of the first of these color structures is given precisely by the set of 'right' diagrams, so it is just 

t One can motivate the argument using string theory, because each primitive amplitude arises from a distinct 
string sector. Assuming that the complete set of string states is inserted between gluon n and fermion 1, then 
A% arises from the Neveu-Schwarz sector and A?- from the Ramond sector. 
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This result holds because when the internal gluon line runs past the n a gluons from SU(N a ), 
and the fermion line runs past the rib gluons from SU(Nb), then the internal gluon must belong to 
SU(N a ), and an extra factor of N a is generated in each such graph. Similarly, in the case that the 
internal gluon line runs past the rib gluons from SU(Nb), an extra factor of Nb is generated in each 
graph, and so the second coefficient is A^. Thus A^ and A^ are separately gauge invariant. This 
argument works even if there are no external gluons between fermions f\ and / 2 , that is, n a = 0. A 
similar argument can be carried out for the fermion (scalar) loop contributions An'^ 1 ^ 2 ^ (respectively 
An'^) simply by adding some extra fermion (scalar) flavors transforming under SU(N a ) but not 
under SU(Nb), in order to separate the L and R diagrams. 



4. Prom Primitive Amplitudes to Partial Amplitudes 

In this section, we present the relation of the qqg . . . g partial amplitudes A n .j to the primitive 
amplitudes An^ J \ An'^ ■ The easiest case is the leading partial amplitude A n; i(lq, 2 q ; 3, . . . ,n), the 
coefficient of the leading color coefficient N c (T a3 . . . T an )^ . Inspecting color flows in the double- 
line formalism, it is easy to see that this amplitude receives contributions only from diagrams 
where the two fermions are adjacent. The A^(l q , 2 q , 3, . . . , n) diagrams contribute with weight 1, 
after noting that the factor of N c in the diagrams supplies the explicit N c in the definition (2.5) 
of Gr^p. The A^(l q ,2 q ,?,^ . . . ,n) diagrams would not contribute at all to A n -\, because of their 
wrong color flow, were it not for the —1/N C term in the SU(N C ) Fierz identity (2.2). Because of 
this term, the A^(l q , 2 q , 3, . . . , n) diagrams contribute with weight —1/N^. The fermion-loop piece 
^j£,[i/2] 2 q ,3, . . . ,n) contributes with weight rif/N c , and similarly for scalar loops. As discussed 

R [1 /2l 

near eqn. (3.4) the primitive amplitude A n (lg,2 g ,3, . . . ,n) vanishes because all its diagrams 
are tadpoles. Putting the contributions together, A n;1 is given by 

'A^iXqi 2q, 3, . . . , n) 

(1.1) 

One can make use of supersymmetry (3.3) to rewrite (4.1) as 
A n;1 (l q , 2,; 3, .... n) = ^ USY (1,-, 2„ 3, . . . , n) - (l + ^\ A*{l q , 2 q , 3, . . . , n) 

+ (l " ^) 41(1* 2 q , 3, . . . , n) + (l + ^ - ^) A° n (l q , 2 q , 3, . . . , n), 

= i 1 + W) A " 2q > 3 ' • • • ' n) ~ JP A ^ SY ^ 2 ~^ • • • ' n ) 

+ ^5) Al(l q , 2 q , 3, . . . , n) + (?LZI* - l^j K (l q , 2 q , 3, . . . , n), 

(4.2) 
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A n -i (l q , 2 q ; 3, . . . , n) — A^(l q ,2 q ,3, . . . ,n) j^A^(l q , 2 q , 3, . . . , n) 

+ ^[^](i_, 2g)3 ,...,n) + ^A^\l qM ...,n). 



where 

Af - -A L 'W - 4 i 't 1 /2] 

n n ' (4.3) 

AS _ /\L,[0] 
n n 

By the expression Af^ SY (lq, 2 q , 3, . . . , n) we simply mean the supersymmetric 5551 ... 5 amplitude 
A^ SY (lg, 2g, 3, . . . , n) with gluino labels replaced by antiquark/quark labels. The combination 

is simpler than J 4^'' 1 ^ 2 ' or An'^°\ because it can be viewed as a chiral matter supermultiplet 
contribution to a supersymmetric amplitude. 

In appendix I we prove that the subleading-color partial amplitudes A n -j > i may be expressed 
as a permutation sum over primitive amplitudes, 



A n]j (l q ,2 q ;3,...,j + l;j + 2,j + 3,...,n) = ^ 

a€COP{a}{(3} 



A L n ^ (<7(l ? ,2„3,...,n)) 



-^^' [1/2] 2„ 3, . . . , n)) - ^A«M (a(l q , 2 q ,3,...,n)) 

(4.4) 

where on € {a} = {j + 1, j, . . . , 4, 3}, ft G {/?} = {1, 2, j + 2, j + 3, . . . , n - 1, n}, and COP{a}{(3} 
is the set of all permutations of {1,2,..., n} with leg 1 held fixed that preserve the cyclic ordering 
of the «j within {a} and of the ft within {/?}, while allowing for all possible relative orderings of 
the ccj with respect to the ft. For example if {a} = {4,3} and {/?} = {1,2,5}, then COP{a}{(3} 
contains the twelve elements 

(1.2.5.4.3) , (1,2,4,5,3), (1,4,2,5,3), (1,2,4,3,5), (1,4,3,2,5), (1,4,2,3,5), 

(4.5) 

(1.2.5.3.4) , (1,2,3,5,4), (1,3,2,5,4), (1,2,3,4,5), (1,3,4,2,5), (1,3,2,4,5) 

(cyclic ordering for a two-element set is meaningless). Note that the ordering of the first set of 
indices is reversed with respect to the second. Formula (4.4) is analogous to the one proven for 
adjoint representation external states in ref. [20]. (In ref. [20] leg n was held fixed; the choice of 
fixed leg is completely arbitrary, and here we find it convenient to hold fixed a fermion leg, labeled 
by 1.) 

The terms independent of nj and n s in formula (4.4) (or equivalently, formula (1.2)) may be 
heuristically understood in terms of 'parent' diagrams, which have no trees attached to the loop, 
as depicted in fig. 5. Performing a color decomposition on ordinary Feynman diagrams, using 
eqns. (2.1) and (2.2), it is easy to see that the set of all parent diagrams feed into both A n .j and 
An'^ in the correct way so that eqn. (1.2) is satisfied for this class of diagrams. Roughly speaking, 
gauge invariance requires all other diagrams to tag along properly with the parent diagrams. The 
right-hand sides of (4.4) and (1.2) do, however, contain diagrams not appearing on the left-hand 
side; in appendix I we prove that all such unwanted diagrams cancel in the permutation sum, and 
also treat the n/- and n s -dependent terms. 

15 



Figure 5. Parent diagrams have no trees attached to the loop. The diagram lines represent 
any particles in the theory. 

Equation (4.4) contains in the permutation sum primitive amplitudes A^ = An'^ with any 
number of gluons sandwiched between the fermions l q and 2 q . If we know the supersymmetric n- 
gluon partial amplitudes, then we can use the reflection symmetry (3.2) plus supersymmetry (3.3) 
to eliminate those A^ with more than (n — 2)/2 gluons sandwiched between l q and 2 q . Below we 
carry out this procedure for the four- and five-point cases. 

As an explicit example, we present the subleading-color four-point qqgg amplitude relations, 
which are particularly simple, 

A 4;3 (lq, 2 q - 3, 4) = A\{\q, 2„ 3, 4) + A^(l q , 2„ 4, 3) + A^(l q , 3, 2 q , 4) + A%{l q , 4, 2 q , 3) 

+ A%(l q , 3, 4, 2 q ) + A%(l q , 4, 3, 2 q ) (4.6) 

= ^ SY (l 5 -,2 5 ,3,4) +Aif Y (l q ,2 q ,4,S) +Aif Y (l q ,3,2 q ,4) , 
where we used eqns. (3.2) and (3.3). In A4.3 the fermion and scalar loop contributions (1.5) vanish. 
The only orderings that do not vanish by equation (3.4) are Af ' [1/2] (l q , 3, 4, 2 q )+Af ,[1/2] (1,, 4, 3, 2 q ), 
and the same combination with [1/2] replaced by [0]; these combinations cancel due to Furry's the- 
orem. In fact the expression for A4.3 reduces to a supersymmetric quantity. One can verify this 
relation using the explicit amplitudes given by KST [16]. 

The five-point relations are of course a bit more complicated, 

A 5 . 3 (l q , 2 q ; 4, 5; 3) = ^ ^{l q , 2 q , (7(3), (7(4), (7(5)) + A%(l q , 4, 5, 2„ 3) + A^(l q , 5, 4, 2 q , 3) 



+ A^(l q , 5, 2 q , 3, 4) + A^(lq,A, 2„ 5, 3) + A 5 L (1 9 -, 4, 2„ 3, 5) + 5, 2„ 4, 3) , 



(4.7) 



4> ;4 (l,-,2,;3,4,5) = 

a€Z 3 
n s ~ Uf 



- A%(l q , 2 q , a(5), a(4), a(3)) - ^(1,-, a(5), 2„ (7(4), (7(3)) 

- (7(5), (7(4), 2„ a(3)) - A^(lq, a(5), (7(4), (7(3), 2,) 
(Al(l q , 2 q , a(3), a(4), (7(5)) + A s 5 (l q , (7(3), 2„ <t(4), (7(5))) 



(4.8) 



+^ 2„ (7(3), a(4), (7(5)) + 4(1,-, (7(3), 2 g , a(4), (7(5)) 

where A§ and A? 5 are defined by equation (4.3). For A 5;3 the n/ jS terms cancel because of equa- 
tion (3.4) plus Furry's theorem; the contribution with ordering (1^,4,5,2^,3) cancels the contribu- 
tion with ordering (l q , 5, 4, 2 q , 3). We have also used the reflection identity (3.2) to convert i?-type 
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fermion/scalar loop contributions into L-type ones. We can further use the reflection identity on the 
J = 1 primitive amplitudes, followed by the supersymmetry identity (3.3) combined with Furry's 
theorem to write these two partial amplitudes as follows, 

A 5;3 (l 9 -,2 9 ;4,5;3) = ^ [A^(l q , 2„ <r(3), a(4), a(5)) + A 5 L (l g -, <x(3), 2„ <r(4), <r(5))] 

- 4 USY (l g -, 3, 2„ 4, 5) - A 5 SUSY (1 ? -, 3, 2„ 5, 4) , 
A 5;4 (l 9 -,2 9 ;3,4,5) = - J] [Af-(1 ? ,2„(7(3),<7(4),<7(5)) +^(1 ?J< 7(3),2„ < 7(4), < 7(5))] 

<tGS 3 



o-eZ 3 



4 USY (l g -, 2„ <7(3), a(4), (7(5)) + 4 USY (1 9 -, ^(3), 2„ a(4), a(5)) 



+ 



+ 1 



(A|(l ff , 2„ a(3), a(4), a(5)) + A?(l ?) a(3), 2 9 , a(4), a(5))) 
(4(1,-, 2„ <7(3), (7(4), a(5)) + 4(1,-, a(3), 2„ (7(4), a(5))) 



(4.9) 

In the next section we give explicit formulae for the primitive amplitudes A SUSY , A L , A s and A? 
with orderings (1,, 2 q , 3, 4, 5) and (1,, 2, 3,, 4, 5). The advantage of the form (4.9) is that all terms 
in the permutation sum may be obtained by a direct relabeling of these primitive amplitudes. 



5. Two-Quark Three-Gluon Primitive Amplitudes 

In this section, we give explicit formulas for all the primitive helicity amplitudes for the two- 
massless-quark three-gluon process qqggg- Using eqns. (4.2) and (4.9), one can form the partial 
color-ordered amplitudes. With these in hand, one can construct either the full amplitude using 
equation (2.4), or the virtual correction to the parton-level cross-section, arising from the color- 
summed interference of the one-loop amplitude with the tree amplitude, using equation (2.11). 

In calculating the primitive amplitudes for A^ ; i(lq , 2+, 3~, 4 + , 5 + ) we used a modification of 
string-based methods for gluons, while for the other helicity configurations we found it convenient 
to calculate in field theory, drawing on lessons from string theory. In particular, we used field-theory 
background-field methods [30] (which are an important ingredient in understanding string-based 
methods in a conventional language [13]). The color-ordered Feynman rules of background- field 
Feynman gauge lead to a number of calculational improvements besides the obvious reduction in the 
number of terms in the three- vertex. In particular, supersymmetric decompositions of the amplitude 
are more evident [3,17,18], and a unitarity-based method [19] that relies on power-counting criteria 
[20,21] is easier to apply. We used unitarity to calculate certain "cut-constructible" contributions to 
the primitive amplitudes A§ (1,, 2, 3,, 4, 5) which enter the subleading-color partial amplitudes (see 
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section 1). To perform the required loop integrations we used the integration methods described 
in refs. [32,33]. 

We present our results in a convention where all momenta are taken to be outgoing, that is 
for the process — > qqggg; helicity conservation along the fermion line thus implies that the two 
fermion legs must have opposite helicity. Our sign conventions for the primitive amplitudes respect 
the antisymmetry of the color-ordered rules in fig. 1 as well as the supersymmetry identities in 
ref. [6]. The overall sign convention for the explicit helicity amplitudes presented here is actually 
opposite to that usually chosen for fundamental representation quarks; however, the overall sign 
of the loop helicity amplitude is irrelevant as long as the tree and loop amplitudes use the same 
convention. An advantage of this choice of signs is that the signs of the tree-level gluino and 
quark partial amplitudes agree, so that the supersymmetry Ward identities hold without any sign 
adjustments. 

We shall express the primitive amplitudes in terms of the Lorentz inner-products Sij = (ki + 
kj) 2 = 2ki • kj, and the spinor inner-products (J I) = = U-(kj)u + (ki) and \j I] = = 

u+(kj)u-(ki), where u±(k) is a massless Weyl spinor with momentum k and chirality ± [8,6]. 

Discrete symmetries reduce the number of independent primitive amplitudes. Parity reverses 
all external helicities in a partial amplitude; it is implemented by the "complex conjugation" op- 
eration "f " , which is the spinor inner-product substitution (j I) — > [/ j] , [j I] — > (I j) , with no 
substitution of i — > —i. For n-gluon amplitudes, parity takes A — > for two-quark (n — 2)-gluon 
amplitudes, with the above sign conventions, there is an extra minus sign, A — > — A^ . Using par- 
ity, we may restrict our attention to two-quark three-gluon amplitudes having either zero or one 
negative-helicity gluon. 

Charge conjugation changes the identity of a quark to an antiquark and vice- versa. Its effect 
on the primitive amplitudes is to exchange the role of 'left' and 'right' contributions, A L <-> A R , 
because the direction of the fermion arrow is reversed. Because the quark and antiquark have 
opposite helicity, charge conjugation allows us to fix the helicity of the antiquark to be negative. 

Finally, equation (4.2) for the leading-color partial amplitude A 5 -i requires color-orderings 
where the antiquark and quark are adjacent, while equation (4.9) for the subleading-color partial 
amplitudes A 5 . 3 and A 5;4 requires also color-orderings with the antiquark and quark separated by 
one gluon. We thus need to present primitive amplitudes for eight distinct helicity /color configura- 
tions: two are infrared-finite, A(— q + q + ++) and A(— q + + q ++); and six are infrared-divergent, 

M~q +q ~ + +), M~q +q + ~ +), M~ q +q + + ~) , M~ q ~ +q + +). M~ q + +9 " +) »nd 

A{— q + + q -\ — ). Only the latter sextet enter into the next-to- leading order corrections to the 
two-quark three-gluon process. 
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We list the required pieces in turn, beginning with the infrared-finite helicity configuration. 
To construct A$.\ we need 



A 



SUSY, 



l g -,2+,3+,4+,5+) = , 



A 5 L (l g -,2+,3 + ,4+,5 + ) 
^(l g -,2+,3 + ,4 + ,5 + ) 



i (12) [2 3] (31) + (14) [4 5] (51) 
32--' 

i 



+ Ag(ir,2+ 3+,4+,5- 



(2 3) (3 4) (4 5) (51) 
(13) [3 4] (4 1) 2 _ (14) (2 4) [4 5] (51) 



[2 3] [2 5] 

48vr 2 I (i 2) (3 4) 2 (4 5) (5 1) ' (1 2) (2 3) (3 4) (4 5) 2 ' [1 2] (3 4) (4 5) 



+ 



+ 



4(lr,2+,3+,4+,5+) = 0. 

For the construction of the subleading-color partial amplitudes A^ and A^ we also need 
4 USY (1^2+,3+,4+,5+) = 0, 

(13) (14) [4 5] (13) 2 [2 3] 



(5.1) 



4(1^2+,3+,4+,5+) = ^-, 



(1 2) (2 3) (3 4) (4 5) (2 3) (3 4) (4 5) (5 1) 



A§(ir,2+ 3+ 4+,5+) = 



+ A§(lr,2+,3+ 4+,5 H 
i (15) (14) [45] 



(5.2) 



, 48vr2 ( i2)(2 3)(4 5) 2 ' 
2+,3+ 4+,5+) = 0. 

The tree amplitude vanishes for this helicity configuration. 

The remaining helicity amplitudes are infrared-divergent, and also require an ultraviolet sub- 
traction. We will present the formulas for unsubtracted amplitudes; to carry out the MS subtraction 
scheme, one should subtract from the leading-color partial amplitudes A^\ the quantity 



cr 



3 1 fll 



2 e 



2 n f 

3iv: 



1 n s 
3N C 



a tree 
^5 



where 



cr 



1 r(l + e)r 2 (l-e) 



(5.3) 



(5.4) 



(4vr) 2 - £ T(l - 2e) 
and D = 4 - 2e. 

We present our results using the dimensional-reduction variant of dimensional regularization, 
with the external gluons treated in four dimensions. This scheme is equivalent at one-loop to 
the string-based 'four-dimensional helicity' scheme of ref. [12]. To convert these results to the 't 
Hooft-Veltman scheme, one must add to A 5; i the quantity 



S 5 = -c r 



1 



1 



\ lIl'IT 



(5.5) 



and modify the coupling constant appropriately [16,12]. We obtained the quantity (5.5) by direct 
calculation in the different schemes, noting that only the singular terms in the integrals contribute 
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to this quantity. This shift is the same as the one found by KST [16] for four-point qqgg amplitudes. 
The universality of the shift S n = — cr |(1 — -jfa ) A^ ee for the two-quark (n — 2)-gluon amplitudes 
A n -i can be inferred from the invariance of physical cross-sections under scheme shifts [16]. It can 
also be inferred from the universality of collinear limits (see eqn. (III.l)). The pure gluon loop 
splitting amplitudes do not depend on which scheme is used [20]; thus in equation (III.l) a shift 
S n -i for °f the form (5.5) implies a shift 5 n of exactly the same form for A n;1 . One may also 

convert the expressions to conventional dimensional regularization. To do so one must account for 
the difference between conventional and 't Hooft-Veltman schemes by having [e]-helicities (gluon 
polarizations pointing into the — 2e dimensions) in observable legs [34]. Since the amplitudes with 
[e]-helicities contain an explicit overall e only the universal poles in e enter and the scheme differences 
may be expected to affect only terms proportional to the tree-level matrix elements. The conversion 
between the various schemes is discussed in ref. [16]. 

The cancellation of infrared (soft and collinear) divergences only occurs after combining the 
virtual corrections presented here with tree-level six-parton contributions to the full next-to-leading 
order process. Various general formalisms exist for constructing infrared-finite distributions numer- 
ically [35,7]. 

For the infrared-divergent amplitudes, it is convenient to decompose the primitive amplitudes 
further in a manner analogous to the decomposition of the five-gluon amplitudes [3] , 

A x = c r (V x Al ree + iF x ) , x = SUSY, L, s, / . (5.6) 

The V factors are purely functions of the momentum invariants s^+i = (ki + ki +1 ) 2 , and do not 
contain other spinor products. All the poles in e are contained in the V factors. There is of course 
some freedom in shifting finite terms between the V and F terms. For the supersymmetric compo- 
nent, the V factor is given by a linear combination of V functions for the all-gluon amplitude [3] 
(after adjusting for the MS subtraction), 

^ SUSY = V&on + ^uon ■ M 

while F SUSY is related to the all-gluon F terms by a supersymmetry Ward identity. The function 
Vj uon is independent of helicities [3] , 

»a~ = 4 £ (-r-) '+i> (-r^) (=^) + !» j ■ (5.8) 

6 ~[ \- S j,j + lJ ~[ \-Sj+l,j+2j \-Sj-2J-lJ 6 

It will be convenient for us to define a related helicity-independent function with the clockwise set 
of double poles from the q to the q omitted, 



9-1 1 / 2 \ e 



(5.9) 
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All the scalar- and fermion-loop primitive amplitudes turn out to be free of poles, so that we may 
take 

V s = V s = . (5.10) 



We denote the dilogarithm [36] by Li 2 , 

Li 2 (x) = 



dy 



ln(l - y) 



(5.11) 



In order to present the results for the remaining functions, we also define the following functions [3], 



Lo(r) 



ln(r) 



Li(r) 



Lo(r) + l 
1-r 



Mr) = 



ln(r) — i(r — 1/r) 
(1-r) 3 



Ls_i(ri,r 2 ) = Li 2 (l - n) + Li 2 (l - r 2 ) + Inn hir 2 - 



7T 



Ls (ri,r 2 ) 
Lsi(ri,r 2 ) 
Ls 2 (n,r 2 ) 
Ls 3 (ri,r 2 ) 



(1 


- n 


-r 2 ) 




l 




(1 


- n 


-r 2 ) 




l 




(1 


- n 


- r 2 ) 




l 




(1 


- ?"i 


- r 2 ) 



[Ls_i(ri,r 2 )] , 

[Ls (ri,r 2 ) +Lo(ri) +Lo(r 2 )] , 



(5.12) 



Ls 1 (r 1 ,r 2 ) + -(Li(ri)+Li(r 2 )) 

Ls 2 (n,r 2 ) + I (L 2 (n) + L 2 (r 2 )) - - -!- 
3 ori br 2 



The Li functions are nonsingular as r — ► 1, and the Ls^ functions are nonsingular as 1 — r\ — r 2 — ► 0. 

For ^4 5 (1^~, 2 + , 3~, 4 + , 5 + ) the functions needed in order to construct the amplitude via equa- 
tions (4.2), (5.6) are: 



.tree = • (1 3)' (2 3) 

5 * (12) (23) (34) (45) (51) ' 



t/SUSY _ yg _ j2_ 
gluon 2 e 



I 1 



-S34 



+ 



2e \-«34 



+ ln 



-S51 

•551 



6, 



(5.13) 



-S12 



3 , 



( 32 ) p/ 

(31) S^°n 



3 (13) (2 3)(41)[2 4] 2L ^i (ggjjgj + 3 (1 3) (2 3) (5 3) [2 5] 2 Ls i ' — — 



(4 5) (51) s 2 51 (34) (45) 



(13)' 



((15) [5 2] (2 3) + (12) [2 4] (4 3)) 



Lo 



-^34 ) 
-«51 / 



'34 



2 (12) (34) (45) (51) V N ' L J x V s 5 i 
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L _ (12) (2 3) (34)(41)[2 4] 3Ls 2 [gf, gg) _ (12) (2 3) (3 5) 2 [2 5] 3 Ls ^ (ggjfg) 
(4 5) (51) ^ (34) (4 5) s| 4 

-■S23 ~^34 \ ,„ „, ,„ „, ,„ , ,9l iC ,, | -Sl2 — S51 

-S34 ' — S34 



(1 3) (2 3) (4 1) [2 4] 2 Ls i (=f£> ^fflj o (1 3) (2 3) (3 5) [2 5] 2 Ls i 
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For use in constructing the subleading-color partial amplitudes, the V and F functions for the 
components of ^(l^, 2~, 3+, 4 + , 5 + ) are 
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The remaining amplitudes are related by discrete symmetries to those presented above; con- 
sider, for example, 
A%(lf, 2* 3" 4", 5"): 

^(1+ 2" 3" 4-,5") =-[Af-(ir,2+ 3+,4+,5 + )] t 

A 5 L (l Q -,2+,3-,4-,5") =A«(2+ 3-,4-,5-,l-) = -^(2+ l-,5-,4-,3-) (5.25) 

= + [^(^,1+5+4+3+)]* • 

The overall signs in these relations drop out in the cross-section, because the signs are the same 
for loop amplitudes and tree amplitudes. 

We have performed a variety of checks on the amplitudes: 

1) a check of collinear factorization for all primitive amplitudes in all channels, illustrated in 
appendix III, providing an extremely stringent check of the primitive amplitudes; 

2) a verification of the supersymmetry identities [15] for the primitive amplitudes ^(l^ , 2^~, 3 + , 4 + , 5" 
and A 5 (l^ , 2+, 3 + , 4~, 5 + ) by explicitly calculating all terms in the supersymmetry relation 
(3-3); 

3) a verification of some of the cuts in amplitudes that were not calculated via cutting methods; 

4) a check on eqn. (4.9) giving A 5 . 3 and A 5;4 , by comparing the poles in e in these quantities 
against explicit formulae for such singular terms in ref. [28]; 

5) a check of formula (2.11) for the virtual part of the color-summed cross-section, by showing 
that the infrared poles in e properly cancel against singular terms in the 2^4 matrix elements 
arising from the integration over soft and collinear phase space [35,7,28]; 

6) a check of the permutation formula for A 5;i in eqn. (4.9), by comparing the fermion loop 
(n/) contributions to the corresponding contribution of the previously calculated [37] process 
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Z — > 37, but with the Z polarization vector replaced with a fermion bilinear (and gluon 
propagator connecting it to the loop); it is not difficult to see that the contributing diagrams 
are identical for these two cases. (The axial coupling does not contribute to Z — > 37 so it does 
not affect the comparison.) We have explicitly verified that the fermion loop contributions on 
the right-hand side of eqn. (4.9) agree with the appropriately modified expressions contained 
in ref. [37] for vanishing fermion masses. 

6. Conclusions 

In this paper, we presented all one-loop QCD amplitudes for two external quarks and three 
external gluons. Combining these results with the ones for five gluons [3] and four quarks and one 
gluon [4], this completes the set of one- loop amplitudes required for calculating next-to-leading 
order corrections to three-jet production at hadron colliders. The computation made use of a num- 
ber of techniques, including spinor helicity [8], color decompositions [5,9,10], string-based meth- 
ods [12,14], supersymmetry methods [15,3,17,18], improved gauge choices [30,31,13], perturbative 
unitarity [20,21], and collinear limits [38,6,23,39]. 

We also introduced 'primitive amplitudes as gauge-invariant building blocks from which ampli- 
tudes containing fundamental representation external legs may be constructed. The usefulness of 
primitive amplitudes follows from their relatively simple analytic structure. In a previous paper 
we obtained a formula [20] valid for adjoint representation states which allows one to obtain all 
subleading-color partial amplitudes from the leading-color partial amplitudes. Using primitive am- 
plitudes we generalized this formula to the case of two external fundamental representation quarks 
and (n— 2)-gluons. Further generalizations to larger numbers of external quarks are straightforward. 

In calculating the amplitudes we made extensive use of supersymmetry identities [15], both 
as a check and to reduce the number of independent amplitudes to be calculated. We verified 
that as the momenta of two adjacent external legs become collinear, the amplitudes presented here 
properly reduce to sums of lower-point amplitudes multiplied by universal splitting amplitudes. 
This provides a stringent consistency check on the amplitudes. 

The one-loop two-quark three-gluon amplitudes presented in this paper constitute, along with 
the five-gluon [3] and four-quark one-gluon [4] ones, one of the major ingredients required for the 
construction of a next-to-leading order program for the prediction of three-jet physics at hadron 
colliders. The infrared singularities in these amplitudes must be cancelled by adding the singular 
contributions from real emission, for example using the formalisms of refs. [7,35]. The program 
also requires the full form of the real emission contributions in non-singular regions, given here 
by the known six-point tree amplitudes [5,6]. Such a program would allow the study of three-jet 
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distributions to next-to-leading order; as with the two-jet case studied extensively by various collider 
detector collaborations, large-statistics data are available. In the three-jet case, one may study a 
richer variety of distributions. The comparison of three-jet rates to two-jet rates, in conjunction 
with such an NLO program, offers the first possibility of measuring the strong coupling constant 
a s in a purely hadronic process deep in the perturbative regime. Beyond three-jet studies, such 
a program also incorporates the elements required to study jet structure beyond the leading non- 
vanishing order available in two-jet NLO programs. 
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Appendix I. Subleading-Color Partial Amplitudes from Primitive Amplitudes 

In this appendix we prove that the subleading partial amplitudes A n -j > \ are given by sums 
over permutations of primitive amplitudes. We will use a similar result for the case of n external 
adjoint particles [20] as an intermediate step in deriving the result for the case that the two external 
fermions are in the fundamental representation. 

First we discuss the distinction between the gauge groups SU (N c ) and U(N C ) = SU (N c ) xU(l). 
If all particles in an amplitude transform as the adjoint representation of SU(N C ), and all vertices 
are given by (combinations of) the structure constants f abc , then there is essentially no distinction, 
because f abc = whenever a corresponds to the £7(1) generator, T au w = 1/y/N^. In other words, 
the U(l) 'photon' is automatically projected out by vertices such as nonabelian vector-boson self- 
interactions. The —1/N C term in the SU(N C ) Fierz identity (2.2) removes the 'photon' explicitly 
by projecting onto traceless hermitian matrices. It can be ignored if only adjoint-representation 
particles are present. 

More generally, the —1/N C Fierz term only affects those diagrams for qqg . . .g where a gauge 
boson propagator is attached at both ends to a line in the fundamental representation — a fermion 
or scalar line. In those diagrams where both of these ends are in the loop, that is with exactly 
one gauge boson propagator in the loop itself (the diagrams contributing to A^(lq, 2 q , 3, . . . , n)), 
the —l/N c term leads to A n;1 contributions only, as discussed above. The only other diagrams 
affected are fermion or scalar loop contributions, where a gauge boson attaches the pinched-off 
external fermion line to the loop. In summary, the nj- and n^-independent parts of the subleading 
qqg . . . g partial amplitudes A n; j > i, and all of the ggg . . . g partial amplitudes, can be analyzed as if 

28 



the gauge group were U(N C ), neglecting the —1/N C Fierz term when working out the color flow in 
the double-line formalism. This result in turn implies that such qqg . . . g double-line diagrams can 
be obtained from a subset of ggg . . . g diagrams by 'color-stripping' — removing a color line that 
flows directly from one gluino to the other, thereby converting the adjoint representation gluino 
into a fundamental representation quark. Consequently, conversion of a subleading-color formula 
for ggg . . . g to the qqg . . . g case is quite straightforward. 

The formula we want to re-derive, and then modify, expresses the subleading-color partial 
amplitudes A^j >2 for ggg . . . g (or ggg ... 5) as a sum over permutations of leading-color partial 
amplitudes A n -i, 



<5(M j l:.y,y • 1 n) = (-I)'" 1 £ <j(a(l,2,...,n)), (1.1) 

a€COP{a}{(3} 



where ai 6 {a} = {j — 1, j — 2, . . . , 2, 1}, £ {(3} = {j, j + 1, . . . , n — 1, n}, and the set of 
permutations COP{a}{/3} is defined below eqn. (4.4). We have added the superscript "adj" to 
(1.1) to avoid confusion with the qqg . . . g amplitudes denoted by A n -j. 

In ref. [20] a string-theory based proof of eqn. (1.1) was presented. We now review this proof, 
but in terms of the color-ordered Feynman rules in fig. 1. The leading-color partial amplitude 
A^{(1, 2, . . . n) associated with Tr(T ai T a2 • • • T a ™) is given by the sum of all color-ordered planar 
diagrams whose external legs follow the cyclic ordering of the color trace, as depicted in fig. 2. 
Examples of five-point diagrams, dressed with color flow lines, are shown in fig. 6. Our convention 
for the direction of the color flow is to follow the reverse ordering of the color trace; this will lead 
to the standard convention for the color arrow following the fermion arrow after conversion to 
fundamental-representation fermions. To compute the subleading-color amplitudes A^j associated 
with the color structure Tr(T ai T a2 ■ ■ ■ T^- 1 )xTr(T^T^+ 1 • • • T a ") directly, one must sum over all 
planar color-ordered diagrams whose corresponding Feynman diagrams can give rise to this color 
structure. These are the diagrams where the cyclic ordering of legs that belong to each trace follows 
the ordering of that trace, but where the ordering in one trace is reversed because the two color 
lines associated with a adjoint particle flow in opposite directions around the loop. For example, in 
fig. 7 the legs follow ordering 123456, but the color structure is Tr(T a4 T a3 T a2 )xTr(T a5 T a6 T ai ). 
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(b) 




(c) (d) 
Figure 6. Examples of the color flow for adjoint representation fermions, using the double- 
line formalism. The Feynman diagrams are gray, and the oriented color lines dressing them are 
black. 



Thus we sum over color-ordered diagrams with the legs permuted over COP{a}{f3}, where 
a = {j — 1, j — 2, . . . , 1}, (3 = {j,j + 1, . . . , n}. We must, however, explicitly exclude one class of 
diagrams whose color flow is incorrect. This is the class of diagrams where indices from both sets 
{a} and {/?} label leaves of the same tree attached to the loop. The color flow in these diagrams 
cannot produce the desired trace structure, because the line attaching the tree to the loop can 
carry only a single pair of color indices, whereas two pairs would be required to join both a and (3 
indices to those elsewhere on the loop. Examples of diagrams contributing to the color structure 
Ti(T ai T a2 )xTr(T a5 T a4 T a3 ) are shown in figs. 6c and 6d, and one should exclude diagrams such 
as the one depicted in fig. 8. (In the string-based derivation [20], another class of diagrams - 
those where a single tree contains all elements of either {a} or {/?} — is also excluded explicitly. 
This class automatically cancels out of the field-theory calculation, so long as the set {a} or {/?} 
in question contains only legs transforming under the adjoint representation.) 
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J 
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Figure 7. The ordering of one trace is reversed as compared to the other. 
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Figure 8. An example of a diagram that does not contribute to the coefficient of the color 
structure Tr (T ai T a2 ) X Tr (T as T a4 T as ) . 

One can divide the set of all diagrams into a 'parent' subset, which have only three-point 
vertices and no non-trivial trees (depicted in fig. 5), and all the remaining diagrams. We refer 
to the latter as 'daughter' diagrams, because each can be derived from some 'parent' diagram via 
a continuous 'pinching' process, in which two lines attached to the loop are brought together to 
a four-point interaction, or further pulled out from the loop, and left as the branches of a tree 
attached to the loop. Repeating this process in all inequivalent ways yields all graphs contributing 
to the same color ordering as the parent diagram, i.e. all of its daughter diagrams. For instance, 
the color-ordered diagram in fig. 8 is a daughter of the parent with ordering 12345 depicted in 
fig. 6a. Daughter diagrams of different parents can be essentially the same diagram. For example, 
if one swaps legs 2 and 3 in fig. 8, one obtains a daughter of the parent with ordering 13245, but 
the two daughter diagrams differ only by an overall minus sign coming from the antisymmetric 
color-ordered three- vertex. Such relations are important for proving eqns. (1.1) and (4.4). 

Let us first focus on the 'parent' subset of diagrams. As mentioned in section 4, by performing 
a color decomposition of ordinary Feynman diagrams and using eqns. (2.1) and (2.2) one can show 
that all 'parent' diagrams feed into both and A^j >1 in the correct way so that eqn. (1.1) is 
satisfied for this class of diagrams. The same arguments apply to those 'daughter' diagrams where 
each pinched-off tree contains only members of the {a} set, or only members of the {/?} set. 

The only thing left to prove is that the class of daughter diagrams specifically excluded from 
A n -,j — where individual trees have both {a} and {/?} members — does yield a vanishing contri- 
bution when summed over the permutations in COP{a}{(3}. Eqn. (1.1) then follows. When such 
diagrams have only three-point vertices on the trees, the diagrams can be arranged so they cancel in 
pairs. The pairs are related by the exchange of an {a} leg with a {(3} leg on a tree; the cancellation 
follows from the anti-symmetry of the three-point color-ordered Feynman vertices in fig. 1 under 
the interchange of the ordering of the two outer legs. For example, the pairs of diagrams in fig. 9 
cancel in the sum. For diagrams with trees containing four-point vertices, using the color-ordered 
rules in fig. 1, the cancellations in the sum over COP{a}{(3} occur in triplets, such as those shown 
in fig. 10. Diagrams with four-point vertices attached to the loop can be decomposed into the 
same color structures encountered above. This provides a purely field-theoretic proof of eqn. (1.1), 
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where all external and internal states are in the adjoint representation, independent of whether 
they are are gluons or adjoint fermions; in particular it applies to the subleading-color ggg . . . g 
super- Yang-Mills partial amplitudes. 




Figure 9. Examples of pairs of diagrams that cancel in the permutation sum. 



Pi 




Figure 10. Diagrams with four-point vertices cancel in the permutation sum in triplets. 

Now consider the modifications necessary for qqg . . . g QCD amplitudes, where the two fermions 
are fundamental representation quarks instead of adjoint representation gluinos. We will show that 
the subleading-color partial amplitudes A n; j > 2(lq, 2 q ; 3, . . . , n) in the full amplitude (2.4), omitting 
for now the n/^-dependent terms from closed fermion or scalar loops, are given by exactly the same 
type of formula as eqn. (1.1), 

A n;j (l 9 -,2 9 ;3,...,j+l;i+2,j+3,...,n)| non/)S = (-l) J " 1 J2 2„ 3, •••,")) , 

a€COP{a}{(3} 

(1.2) 

where on G {a} = {j + 1, j, . . . ,4, 3}, ft G {(3} = {1,2, j + 2, j + 3, . . . ,n - l,n}, COP{a}{(3} is 
defined in exactly the same way as for the adjoint representation formula, and ny s means either 
n s or rif. 

The main difference between the adjoint formula (1.1) and the fundamental formula (1.2) is 
that the adjoint formula lumps 'left' diagrams, where the fermion goes around the loop on the left 
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side, together with otherwise identical 'right' diagrams, whereas the fundamental formula keeps 
the two separate, as required by their different color flows. To derive (1.2) from (1.1) we remove 
a single color line from a special subset of the ggg . . . g partial amplitudes A^- on the left-hand 
side of (1.1), those where the two gluino charge matrices are in the same trace and are adjacent to 
each other; this subset is in one-to-one correspondence with the qqg . . .g partial amplitudes A n -j. 
We show that this color-line removal corresponds to dropping the unwanted 'right' set of diagrams 
on the right-hand side of (1.1), thus converting A^{ to A^. As discussed above, the —l/N c Fierz 
corrections can be ignored here. 

Let us focus on the coefficient of the same-trace color structure Gr^ ;j - (for some fixed j) in the 
modified color decomposition (2.8), 

Tr(T a3 . . . T a ' +1 ) Tr(T fl T h T a i+ 2 ...T a "), j = 3, ...,n-2, 

(1.3) 

Tr(r 3 ...T a ") Tv(T fl T h ), j = n-l, 

whose corresponding partial amplitude is A*.j(3, . . . 1§, 2 g , j + 2, j + 3, . . . , re). Our convention 

for drawing the color-dressed diagrams is to impose a clockwise ordering on the legs associated with 
the color trace containing the fermion color matrices and a counterclockwise ordering on the legs 
associated with the other trace. In every double-line diagram contributing to this partial amplitude, 
there is a color line that runs directly from gluino 2 to gluino 1 along the gluino line, as examples 
in fig. 6 illustrate. If we remove this color line (as shown in fig. II), then (1.3) is converted to 

Grg ) (3,...,j + l;j + 2,...,n) = Tt(T a ^...T a ^) (T**+» . . .T°«)£ , j = 3,...,n-2, 

0^(3,..., n) = Tr(T^...T a ") 5£ , 

(1.4) 

that is to the fundamental representation color structures given in eqn. (2.5), with coefficients A n .j. 




(b) 




(c) 



(d) 



Figure 11. Examples of color stripping; the color lines running directly from gluino 2 to 
gluino 1 in fig. 6 have been removed. 
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A given ggg . . . g diagram, with one color line running between the two gluinos removed, can be 
interpreted as a qqg . . . g diagram, but only if the color line to be removed runs along the fermion side 
of the loop. This requirement eliminates half the diagrams, namely the 'right' diagrams contributing 
to A]%, leaving precisely the desired 'left' diagrams contributing to A%. As a particular example, 
the two routings of the gluino through the diagrams depicted in fig. 12a, b contribute to the color 
structure (1.3). But in fig. 12b the color line running directly between the two gluinos must run 
along the gluon side of the loop in order to generate (1.3); therefore this diagram should be dropped 
in converting to qqg . . . g, while fig. 12a should be kept. This shows that replacing A n .\ by A^ on 
the right-hand side of eqn. (1.1) is the correct prescription for the 'parent' diagrams. However, we 
must again show that all unwanted diagrams with attached trees cancel in the permutation sum. 
As was the case for the adjoint representation case, in each term on the right-hand side of eqn. (1.2) 
we are including diagrams which do not belong on the left-hand side, because the color flow they 
represent does not allow them to contribute. The argument is similar to the one for the adjoint 
case. If diagrams contain trees with leaves labeled by indices from both sets {a} and {/?}, then 
such diagrams cancel in the permutation sum exactly as for the adjoint case. These cancellations 
are due to the antisymmetry of the color-stripped vertices, and are thus independent of the color 
representation of the fermions. Note that 'left' diagrams cancel against 'left', and 'right' against 
'right'. This completes the conversion of the adjoint formula (1.1) to the fundamental formula (1.2), 
excluding the contributions of closed fermion or scalar loops. 




Figure 12. Both 'left' (a) and 'right' (b) diagrams contribute to adjoint representation 
gluino partial amplitudes, but only the 'left' diagram contributes for fundamental representation 
quarks. 

We turn next to the contributions of closed fermion or scalar loops in the fundamental repre- 
sentation. In this case, neglecting the —1/N C term in the SU(N C ) Fierz identity (2.2) leads only 
to the n/ jS -dependent terms in A n -i in eqn. (4.1), and does not yield a contribution to A n .j > i. 
Including the —1/N C , or U{\) subtraction, term in the gluon propagator connecting the external 
qq line to the fermion/scalar loop decouples the color flow for the tree containing qq from the loop 
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color flow. However, the color-stripping argument can still be applied to these terms, if we allow 
the stripped color line to propagate down the 'photon', around the loop, and back again along the 
'photon'. We start with the ggg ■ ■ -g double-line configuration where a color line starts at gluino 
2, flows down the gauge boson line connecting the external fermion legs to an adjoint fermion (or 
scalar) loop, flows around the loop and then returns through the gauge boson line to wind up at 
gluino 1, as shown in fig. 13a. We remove this color line, as shown in fig. 13b, to obtain the 
U(l) subtraction contribution for qqg...g. In this case it is the 'right' type diagrams that have 
the correct color flow to be strippable, so the appropriate formula (including the Fierz factor of 
-l/N c ) is 

An;j (lgj 2q,; 3, . . . , j + 1; J + 2, J + 3, . . . ,n)\ nf3 terms 



(-ir 1 E 

a£COP{a}{(3} 



{a(h: 2g , 3, . . . , n)) - ±A*>W ^(1,, 2 g , 3, . . . , n)) 



(1.5) 

The cancellation of unwanted diagrams with attached trees from the right-hand side of eqn. (1.5) 
works as in the no-ny )S case. 





(a) (b) 

Figure 13. Color stripping the gluon line connecting the external fermion line to the fermion 
loop in (a) produces the 'photon' subtraction diagram in (b). 



Thus the final formula for subleading-color qqg . . . g partial amplitudes in terms of primitive 
amplitudes is the one given in eqn. (4.4). 



Appendix II. Four-Point Amplitudes 

In this appendix, we collect the four-gluon and two-quark two- gluon amplitudes, needed for 
checking the collinear limits of the two-quark three- gluon amplitudes. These amplitudes agree with 
the results of KST [16]. Note, however, that these authors used a different color decomposition and 
overall sign convention than used in this paper. 

We begin by listing all tree amplitudes that appear in the collinear limits, eqn. (III. 1) , of the 
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five-point amplitudes presented in this paper, 

< ee (l-,2-,3 + ,4+) = i 
< ee (l-,2+,3-,4+) = i 
4 ree (l g -,2+,3-,4+) = i 
4 ree (l g -,2+,3+,4-) 
< ee (l g -,2-,3+,4+) 



(12> Z 



(12) (2 3) (34) (41)' 



(IS)* 



(12) (2 3) (34) (41)' 



(1 3y (2 3) 



= i 



(12) (2 3) (3 4) (41) 

(14) 3 (24) 
(12) (2 3) (3 4) (41) 

(12) 3 (3 2) 



(HI) 



(12) (2 3) (3 4) (41) 

Now consider the one-loop four-point amplitudes, beginning with the four-gluon amplitudes. 
Amplitudes with only external gluons may be decomposed in terms of contributions to supersym- 
metric multiplets 



^<l.V..,») = *+(4-£)4 + (l + £-£)*, 



(11.2) 



where A 9 is the contribution of an N = 4 multiplet, — A? the contribution of an N = 1 chiral 
multiplet, and A s the contribution of a complex scalar. 
The finite four-gluon amplitudes are 



Ai(l+ 


,2+ 


3+ 


4 + ) 




,2 + 


3+ 


4 + ) 


Al(l+ 


,2 + 


3+ 


4 + ) 


4(1" 


,2 + 


3+ 


4 + ) 


4(i- 


,2 + 


3+ 


4+) 


Al(l~ 


,2 + 


3+ 


4 + ) 



•S12S23 



48vr 2 (12) (2 3) (34) (41) 



(II.3) 



0. 



(2 4) [2 4]' 



48tt 2 [12] (2 3) (34) [41] ' 

The tree amplitudes vanish for these helicity configurations. 

For the amplitudes containing infrared and ultraviolet singularities, we further decompose the 
amplitude into V and F pieces of the sort given in eqn. (5.6) 



Al = c r (V x A¥* + iF x ), 



x = g,f,s. 



The universal functions 



V 9 

4,gluon 



-S12 



+ 



-S23 



+ In 



'{—) 



(II.4) 



(II.5) 
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appear in the four-point amplitudes with infrared divergences. 
For the four-gluon amplitude At ; i(l~, 2~ , 3 + , 4+) we have 

F 9 = 0, 

I ( ^ 



V 9 = V 9 

4,gluon ' 



V s 



-S23 

V s 2 
1~ + 9' 



F f = Q, 
F s = 0. 



For A 4; i(l-,2+,3-,4+) we have 



= v? , 

4,gluon ' 



1 



/i 



2 \ e 



2e VV- s i2 
V' 2 
1" + 9' 
<13> 4 



+ 



/i 



-S23 



F 3 = , 

2 \ e \ 

2. 



(12) (2 3) (3 4) (41) 



(13)' 



2s 2 

S12S23 



fln 2 f^Uvr 2 U- gl2 ~ g23 In'"' 12 

V V-S 23 ^ / 2 S13 

(— ) 



(12) (23) (34) (41) V s 2 13 



+ 



-S23 
S12S23 



-•S23 



'13 



gig ~ ^23 
S13 



1 + 



'13 



6S12S23 



In 



+ 1 



For the finite amplitudes with two quarks and two gluons we have 
^ SUSY (ir,2+,3+,4+) = 0, 

i (12) [2 4] 



Af(ir,2+,3+,4+) = - 
AS(ir,2+ 3+,4+) = - 



32vr 2 (2 3) (34) 

i s 23 (12) [2 4] 
48vr 2 si 2 (2 3) (3 4) 



+ ^(ir,2+ 3+,4+), 



i (13) [2 4] 
32vr 2 (2 3) (34) 



(II.6) 



(II.7) 



(II.8) 



A{(lr,2+ 3+ 4+)=0. 

Af SY (lr,2+,3+,4+) = 0, 
A 4 L (l g -,2+,3+,4+) 

A2(ir,2+,3+ 4+) = 0, 

A{(ir,2+ 3+ 4+) = 0. 
For the helicities containing infrared and ultraviolet singularities we decompose the amplitude 
as in the five-point case (eqn. (5.6)), and again we have V s = V- f = 0. For A 4 (l^, 2+, 3~ , 4 + ) we 



(II.9) 



have 



ySUSY = yg _ 

4,gluon 2g 



-S12 



+ 



-S23 



-6. 



+ ^7 



1 / ^ 2 



4,gluon ' £ 2 ,_ Sl2 



2e V-S12 



5 

2 ' 



(11.10) 
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^SUSY 



= 3 



(1 3} (2 3) 



S12S23 



- 2 (5S)- 2 ) 



1 S12 - S23 



F L = 



(12) (2 3) (34) (41) [ 2s 2 13 V" \-s 23 J ' " J + 2 s 13 ln 
(13) 3 (2 3) S12 



(— ) 



(1 2) (2 3) (3 4) (41) 2s 13 
F s = 0, 

= 0. 



For At(ir,2+,3+,4-) we have 



«13 



-S23 



s 13 S 13 



-«23 



+ IT 2 + 1 



(11.11) 



t/SUSY _ 

4,gluon 



-S12 



6, 



"4,gluonT £ 2 I 



A 4 



2e V-S12 



(11.12) 



nSUSY 



= 0, 



F 



(14) 3 (2 4) gig 
(12) (2 3) (3 4) (41) 2s 13 
F s = 0, 

= . 



In" 



(— ) 



(11.13) 



For At(l^,2-,3+,4+) we have 



t^SUSY = v g 

v M,gluon 



yL _ _ySUSY 
2 



e \-S23 



6, 



(11.14) 



F SUSY = 0, 



F L = 0, 



F s = , F f = 0. 

The remaining amplitude, A^l^" , 2+, 3+,4~), may be obtained via a reflection, 



(11.15) 



A 4 (l g -,2 + ,3+,4-)=^4(l Q -,4-,3+,2+) 



(11.16) 



Appendix III. Collinear Limit Checks 

In this appendix, we illustrate the use of collinear limits in verifying the correctness of explicitly 
calculated amplitudes, including signs and normalizations. As the momenta of two color- adjacent 
external legs become collinear, the amplitudes factorize into sums of lower point amplitudes mul- 
tiplied by 'splitting amplitudes'. The constraints provided by the collinear limits are sufficiently 
restrictive that they can be used to construct ansatze for higher-point amplitudes based on known 
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lower point amplitudes [23,20]. The collinear-limit checks may also be used in numerical programs, 
when converting the amplitudes presented in this paper to physical cross-sections. 

At one loop, the collinear limits of color-ordered one-loop QCD amplitudes have the form 



where k a — > zP and kb — > (1 — z)P, with P = k a + kb, P 2 = s a b — ► 0. The tree and loop splitting 
amplitudes, Split ^ e and Split '°^ p , behave as 1/y/Sab in this limit. This formula holds for any 
of the primitive amplitudes — which at loop level may carry the additional labels J = 1,1/2,0 
(for n-gluon amplitudes), or x = SUSY, L, s,f (for two-quark (n — 2)-gluon amplitudes) — as 
well as for the leading-color amplitudes A n -i, provided that one assigns the correct labels to all 
'loop' quantities in (III. 1) . A tabulation of the splitting amplitudes appearing in massless QCD 
computations was given in appendix B of ref. [20]. The results were given in terms of leading-color 
partial amplitudes A n -i rather than primitive amplitudes; below we shall convert the results to the 
latter form. We describe the collinear behavior of amplitudes before subtraction of the ultraviolet 
pole. A proof of the universality of the splitting amplitudes, limited to external gluons and scalar 
loops, was outlined in refs. [23,39]; a more complete treatment will be given elsewhere. The power 
of the collinear-limit constraint arises from the fact that relationship (III. 1) must hold in every 
channel. 

In order to use the splitting amplitudes in appendix B of ref. [20] for the primitive amplitudes 
presented here, we must take into account that the overall sign convention used in ref. [20] for 
the helicity amplitudes A(lq , 2+, 3, . . . , n) is opposite to that used for primitive amplitudes in the 
present paper. (The overall sign convention for A(l^~,2~,3, ... ,n) is the same.) With the former 
sign convention, one cannot relabel q — > q, q — > q in amplitudes without introducing extra signs. 
Equivalently, one cannot directly interpret the tree-level quark amplitudes as gluino amplitudes, 
because the extra constraints on gluino amplitudes imposed by the Majorana nature of the gluino 
are not satisfied with the former choice of sign. The change in amplitude conventions implies a sign 
change in the second and fourth tree splitting amplitudes in eqn. (B.5) of ref. [20]. With these sign 
changes the tree splitting amplitudes in that appendix may be used for the primitive amplitudes 
presented in this paper; furthermore, one may relabel q — > q, q — > q in the splitting amplitudes to 
get the additional orderings Split_ A (g, g) and Split_ A (g, q) which appear in limits of the primitive 
amplitudes. 



[23,20] 




(III.l) 
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For all splitting amplitudes which do not vanish at tree level, the proportionality constant r$ 
is defined by 



Split 1 °° p (a A %6 Ai ') = c r x Split*if (a A - , ^ ) x r s {-X,a K ,b Xb ) . (III.2) 

For any collinear limit of a supersymmetric primitive amplitude A^ USY , rs is given by the function 
r| USY given in eqn. (B.13) of ref. [20]. 

For the non-supersymmetric primitive amplitudes (x = L, s, /) we first discuss the case g — ► gg. 
Because the fermion line is routed on a definite side of the loop in these components, the loop 
splitting amplitudes depend on whether the pair of adjacent collinear gluons a, b are between q and 
q (A*(lq, . . . ,a, b, . . . , 2 q , . . . , n)) or between q and q (A*(lq, . . . , 2 q , . . . , a, b, . . . , n)). In the first 
case the x = L,s,f loop splitting amplitudes all vanish. In the second case they are given by the 
pure-gluon (J = 1) contribution for x = L, the negative of the sum of the J = 1/2 and J = 
contributions for x = /, and the J = contribution for x = s: 

( 0, x = SUSY, / ; 

Splits (a+ 6+) = ' (III.3) 

{ Split^(a + ,5 + ), x = L,s 

(this is the one case where the tree splitting amplitude vanishes), and 

r L s {±,a,b) = rj! ] (±,a,6) , 

r f s (±,a,b) = 0, (III.4) 
r> s (±,a,b) = rj° ] (±,a,&) . 

The quantities Split^' (a + , b + ) and r| 7 ' are given in eqns. (B.8) and (B.9) of ref. [20]. 
The rs functions for the non-supersymmetric g — > qq splitting amplitudes are 



r L s (±,q^,q ± ) = -1 



9 \ e / 9 

A* \ ( V 



r s s (±,<l T ,q z 



e 2 [\z(l - Z)(-Sq q )j \(-Sgq) 

13 / fi 2 \ e , . , , . , vr 2 83 
+ 6^((4-)) +21n(.) ln(l -,)-- + - 

^ 2 



(III.5) 



3c \ ( Sqq 



For the non-supersymmetric q ^ qg splitting amplitudes, one obtains different results de- 
pending on whether the gluon (denoted by a) is before or after the quark, with respect to the 
clockwise ordering of the vertices. The difference is again due to the routing of the fermion line in 
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the primitive amplitudes. We have 

r^(q-,a+) = rf(g + ,a") = f L (l-z,s qa ), 
r L s {q-,a-) = r^(q+,a+) = f L (l-z,s qa ) + 



(III.6) 



r%(a + ,q+) = r%{aT ,q~) = f L {^s aq ) + -, 

r s( a ~^ + ) = r s( a+ >9") = f L (z,S a q) , 

rs(a + ,q-) = r^{a~ ,q+) = f R (z,s aq ), 
r L s (a-,q-) = r^{a+,q+) = f R (z,s aq )-^, 

r L s {q + ,a + ) = r L s {q~ ,a~) = f R (l- z ,s qa )-^^, 
rl(q + ,a-) = r L s {q~ ,a+) = f R (l - z, s qa ) , 

where the functions f L and f R correspond to the leading- and subleading-color parts of the function 
/ defined in eqn. (B.ll) of ref. [20], 

/ L ( z .») = -3(^5)' -"»(»-*)• 

/*<*•> = - Miri))' + M(5))'- Li2(3) ' (nL7) 

f(z,s) = f-(z, ll )--Lf*(z,s). 

C 

The x = s, f parts of the (/ — > qg splitting amplitudes all vanish. 

As an illustration of the collinear limits consider the amplitude A[ 2+, 3~, 4 + , 5 + ) given 

by 



AH\- 2+ 3" 4+ 5+) = icrF' - ia (13) t 45 ] V-^/ (IIL8) 
A 5 (l,,2 g ,3 ,4 ,5 j _ ic v b - *c r (i2)(45) ^ 4 _ si2 . 

The collinear limits of may be considered independently of all other primitive amplitudes because 
they are separately gauge invariant; in QCD amplitudes A J 5 enters with a coefficient proportional to 
the number of fermions rif, which distinguishes it from the other primitive amplitudes. (Although 
Al also enters QCD amplitudes with coefficients containing n/, it does so only in the combination 
n s — rif and it may therefore be treated independently.) 

First consider the limit as gluon 3 becomes collinear with gluon 4 (3 || 4), and also the 
4 || 5 collinear limit. In either of these cases, eqn. (III. 8) does not contain a collinear singularity. 
Compare this result to the expectation from eqn. (III.l). These limits are particularly simple 
to analyze because the A{ amplitudes with two quarks and two gluons vanish for all helicities, 
A^(l q , 2 q , 3,4) =0 (see eqns. (II. 8), (11.11) and (11.13)). Additionally, the loop splitting amplitudes 
vanish for all helicities, Split£(a ± , 6 ± ) = 0, where Split£ = -(Split^ /2] + Split^ ] ), and Split^ /2] , 

41 



Splitj^ are defined in appendix B of ref. [20]. Thus, all contributions to the right-hand side of 
eqn. (III. 1) vanish for the 3 || 4 and 4 || 5 collinear limits, in agreement with the lack of collinear 
singularities in these channels in the amplitude (III. 8). 

The lq || 2 q collinear limit is a bit less trivial since the expression (III. 8) does contain a collinear 
pole. Using the explicit value of A[ in eqn. (III. 8) we have 



i||2 . z (P3) 2 [4 5] ln (=fi) 



^ -^ "; 5 f (....9) 

where we have taken k\ = zP and k<i = (1 — z)P. Now compare this to the result obtained from 
the collinear formula (III.l). The necessary four-gluon amplitudes are 



^ ree (P + ,3-,4+,5 + ) = 0, 
< ee (F-,3-,4+,5+) 



(P3) 



4 



(111.10) 



(P3) (34) (45) (5P)' 
4(P+,3-,4+,5+) = 0, 

4(F-,3-,4+,5+) = -c r Q +2^ < ee ( J P-,3-,4+,5+), 

while the relevant splitting functions are 

S P litrr,g+) = ^, 

/1 / 2 \e \ (111.11) 

s P .,^-., + ,^ r (-(^-) + S )JL. 

Plugging these results into eqn. (III.l) we obtain for the non- vanishing terms 
4(ir,2+,3-,4+ 5+) 

^ Splitf e (lr,2+) 4(p-,3-,4+,5 + ) + Spli<(ir,2+) < ee (p-,3-,4 + ,5 + ) (nL12) 

%CT (12) n U S34 y'(P3)(34)(45)(5P) 

which reproduces eqn. (III. 9), after spinor helicity simplifications. 

One can continue in this way, systematically verifying that the collinear limits in all channels 
are correct. We have done so for all amplitudes presented in this paper. 

In checking the collinear limits, the behavior of most of the functions appearing in the ampli- 
tudes is straightforward to obtain. We present here only the one function with a slightly complicated 
limit. We take the collinear limit of two color-adjacent momenta k c and k c+ i; denoting the sum by 
P, the momentum fraction z satisfies k c = zP and k c+ \ = (1 — z)P in the limit. Furthermore, we 
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will denote the limit of s c+ 2,c-2 by s and the limit of s c _2, c -i by t. With these definitions we have 

T / — Sc-2,c-l — «e + 2,c-2 \ 

2 V ' -»c,c+l J C||C+1 1 1 



6 c,c+l 



2st s 
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>c,c+l 



+ 



-t / Lt(s + t) 2 2t 2 (s + t)_ 
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,s(s + t) 2 2s 2 (s + t)_ 



1 1 

+ 



2t 2 (s + t) 2s 2 (s + t) 2st(s + t) ' 



(111.13) 



Appendix IV. Mixed Photon-Gluon Amplitudes 

The same primitive amplitudes used to construct the two-quark (n — 2)-gluon amplitudes can 
also be used to construct amplitudes with one to (n — 2) photons replacing gluons. The two-quark 
(n — 2)-gluon color decompositions (2.3) and (2.4) are valid for external gauge bosons in U(N C ) 
as well as SU(N C ), because the Fierz subtraction term in the gluon propagator (eqn. (2.2)) does 
not contribute unless the gluon is sandwiched between two fundamental representation lines. Thus 
partial amplitudes with photons may be obtained by substituting the photon generator matrix 
(which is proportional to the identity matrix) into the appropriate color decomposition formula, 
and grouping together terms with the same color structure. We illustrate the construction explicitly 
for two-quark one-photon (n — 3)-gluon amplitudes. 

At tree level, amplitudes with one photon have a color decomposition similar to that of pure 
nonabelian ones, 

^i 7 tree = Q ^2 e5 n- 3 £ (T°'(») ■ ■ ■ T^n-i))^ «A^ tree (l ff , 2,; (7(3), . . . ,<7(n - l);n) , (IV.l) 

where e is the QED coupling constant, and Q is the charge of the quark. The photon is taken to 
be the last leg, n. The partial amplitude Aj^ tree is related by a 'decoupling' equation [6] to the 
pure nonabelian partial amplitudes A^ ee , 

Al? tree (l q , 2,; 3, . . . , n - 1; n) = < ee (l g -, 2 q ; n, 3, 4, . . . , n - 1) + < ee (l,-, 2 q ; 3, n, 4, . . . , n - 1) 

+ --- + A t n ree (l,-,2,;3,4,...,n-l,n) . 

(IV.2) 

Equation (IV.2) is obtained by substituting the photon generator matrix T a ™ oc 1 into eqn. (2.3) 
and collecting terms. Note that unwanted diagrams on the right-hand side of (IV.2) — those 
coupling the photon to gluon lines — cancel out in the sum. 
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We can perform similar decompositions at one loop. We keep only 0(e) contributions; we 
do not include the C(e 3 ) contributions where an internal photon line is exchanged between two 
charged lines. If we set aside the n/- and n s -dependent pieces, then the color factors are again 
simply the color factors for the two-quark (n — 3)-gluon amplitude (2.5), 



n-2 

.(99) 

v / V ... ! < / / / V T.I 

\n f 

i=i aes„- 3 /s„-i: 



<% ,n s =o = Q^9 n - X £ £ Gr H>( 3 ' • • • > «-!)) ^-(Iff. 2 95 *(3, • • • , n - 1); n) 



,n s =0 

(IV.3) 

and the partial amplitudes with one photon are again given by sums over the two-quark rest-gluon 
partial amplitudes, inserting the photon in all inequivalent gluon locations, 

2 q ; 3, . . . , j + l;j + 2, . . . , n - 1; n) 

71 f ,7i s =0 

= A n . j+1 (l q ,2 q ;n,3, . . . , j + 1; j + 2, ... ,n- 1) H h 2 9 ; 3, . . . ,n,j + 1; j + 2, . . . ,n - 1) 

+ ^ n;i (lg,2 9 ;3, . . . , j + l;n, j + 2, . . . ,n - 1) H h 2 q ; 3, . . . , j + 1; j + 2, . . . , n - 1, n) ; 

(IV.4) 

the unwanted diagrams again cancel. 

For those contributions that do contain a (charged) fermion or scalar loop, we must in effect 
consider separately the diagrams where the photon couples to the external fermion line, and the 
diagrams where it couples to the closed fermion (or scalar) loop, since these different fermion lines 
may have different charges. This can be accomplished (for j > 1) by considering separately the 
contributions where the photon replaces a gluon within the first set, 3 . . . j + 1, associated with the 
trace of gluon matrices in Gr^{ •, or within the second set, j + 2 . . . n, associated with the string of 
gluon matrices whose i2,H component appears in the color factor. In the former case, the factor of 
Q along with rif (or n s ) will be replaced by the trace over the fermion charge matrix (respectively 
the scalar charge matrix); in the latter case, the amplitude will continue to appear with factors of 
Q and rif (or n s for scalar- loop contributions). 

More concretely, define A^/f^ to be the coefficient of nf/N c in A n .j, 

A^!J ] 2 q ; 3, . . . , n - 1, n) = N c —A n , 3 {l- q , 2 q ; 3, . . . , n - 1, n) . (IV.5) 

Using equations (4.2), (4.4), we may relate A^/^ to primitive amplitudes, 

4!/i 2] 2 g ; 3, . . . , n - 1, n) = -A f n (l q , 2 q ; 3, . . . , n - 1, n) - A s n (l q , 2 q ; 3, . . . , n - 1, n), 
4 1 i g 1 (l ff> 2,;3 > ...,n-l > n) = (-1)' £ A^ (a(l q , 2 q , 3, . . . , n)) . ( IV - 6 ) 

aeCOP{a}{/3} 

We can then construct the two different parts of the one-photon partial amplitude for each j > 1, 
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as described above, 

A%f /2]:I 0-v 2,; 3, . . . , j + 1; j + 2, . . . ,n - l;n) 

= A^m^nA . . . , j + l;j + 2,... ,n- 1) + • ■ ■ + ^.^(l,-, 2 9 ; 3, . . . ,n, j + 1; j + 2, . . . ,n - 1) , 
4y- [1/2]: "(l?> 2 g ; 3, . . . ,j + 1; j + 2, . . . , n - 1; n) 

= A^.f (1?,2,;3, . . . , j + l;n,i + 2, . . . ,n- 1) + ■ ■ ■ + A^f 2 g ; 3, . . . , j + 1; j + 2, . . . , n - 1, n) . 

(IV.7) 

For j = 1 the color flow is different, and we get instead 

KT^' 1 ^ 2 9 ; 3, . . . , n - 1; n) = ^.f (1,-, 2,; n, 3, . . . , n - 1) + ■ ■ ■ + <( 2] (1,-, 2,; 3, .... n - 1, n) 

+ ^[ 1 /2] (1 _ ;n)2(j)3) ... )n _ 1)) 

^ [1/2]:/7 (l„2 g ;3,...,n-l;n) = -A^ 1 / 2 ^-, n, 2„ 3, . . . , n - 1). 

(IV.8) 

Using these pieces, we write out the ^/-dependent pieces in the full amplitude, 
Al% = V2eg^ g E GrJS>(3, . . . ,n - 1)) ^^'(l^aft . . . ,n - l);n) 

i=l c6S„_ 3 /S„_i ; ,,- 

n n-2 

+ ^ v^" 1 E E Gr -1' W 3 ' !)) ^ [1/2]:J7 (1^ 2 9 ; <t(3, . . . , n - 1); n) , 

i=i ires„_3/s„_i ; j 

(IV.9) 

where Q/ is the fermion charge matrix, and Tif represents the trace over flavors. An analogous 
decomposition holds for the contributions proportional to n s . 
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